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The theory of the hydraulic analogy, that ic, the analogy 
between vat or flow with a free surface and two -dimensional compressible 
gas flow, and the limitations and conditions of the analogy are 
discussed. A test was run using the hydraulic analogy as applied 
to the flow about circular cylinders of various diameters at subsonic 
velocities extending into the supercritical range. The apparatus and 
techniques used in this application are described and criticized. 
Reasonably satisfactory agreement of pressure distributions and 
flow fields existed between water and air flow about corresponding 
bodies. This agreement indicated the possibility of extending 
experimental compressibility research by new methods. 


INTRODUCTION 


An analogy exists between waller flow with a free surface end 
two-dimensional compressible gac flow (hydraulic analogy) . The 
water- must flow over a smooth horizontal surface bounded by 
vertical walls geometrically similar to the walls bounding the 
corresponding compressible gas flow. 

The mathematical basis of this hydraulic analogy was presented 
by Riabouchinsky in reference 1, in which lie also described his 
apparatus for investigating the flow in a Laval nozzle. In 
reference 2, he extended the theory .to include drag considerations 
and outlined the probable usefulness of the hydraulic analogy'. 
Binnie and Hooker, in reference 3> obtained surveys along the 
center line of a channel with a constriction. By employing the 
characteristics method to calculate accurately the flow in a Laval 
nozzle, Preiswerk, in reference 4, demonstrated conclusively that 
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the methods of gas dynamics can be applied to water flow with a 
free surface. 

The National Advisory Committee for Aeronautics became interested 
in the hydraulic analogy because it seemed an easy and inexpensive 
way of studying two-dimensional compressible gas flow; in particular, 
phenomena occurring in air at speeds too high for visual observations 
could be observed at very low speeds (3 or k fps) in a water channel- 
Preliminary investigations were made in the Langley tank no- 1, 
where difficulty with the vertical accelerations - assumed negligibly 
small in the analogy - was experienced. 

A water channel was designed and constructed in the Langley 8-foot 
high-speed tunnel in the spring of 19*4-0. The channel was so 
constructed that flow fields involving both subsonic and supersonic 
velocities about aerodynamic bodies could be investigated. The 
value of the analogy in such flow fields has not been previously 
demonstrated. The development of the measuring apparatus and 
techniques is presented herein. The application of the analogy to 
flows through nozzles and about circular cylinders at subsonic 
velocities extending into the supercritical range is also presented. 


SYMBOLS 


c v 

7 

T 

P 

p 

h 

a 

a 

P 

h 


specific heat at constant pressure 

specific heat at constant volume 

adiabatic gas constant, ratio of c. 0 to c y 

absolute temperature of gas 

mass density of gas 

pressure of gas 

enthalpy or total heat content 
dynamic pressure of gas (jjpV 2 y 
surface tension of liquid 
pressure coefficient 
viscosity of liquid 



of gas 
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a. speed of sound in gas 

V velocity of flow 

d water depth 

M Mach number, atresia value unless otherwise indicated 
( — • for gas; — for wuteA 

V .& VC'd j 

g acceleration of gravity 

F c compressibility factor 

x, y rectangular coordinate cocos 

u, v components of velocity in x and y directions, respectively 
^ velocity potential in two-dimensional flow 



X 

U 


V 


B 


A 


wave length of surface waves in fluid 

velocity of propagation of surface waves in fluid 

diameter of circular cylinder 
/VLp \ 

.Reynolds number l j 

* p* 

angular measurement clockwise about cylinder; at stagnation 
point, 0=0° 


Subscripts: 

No subscript any value of variable 

0 value at stagnation (V = 0) 

st value for completely undisturbed static conditions, 

no flow in channel 

1 local valuo of variable; value at surface of model, 

at channel walls, br in field of flow 

s value in undisturbed stream 


1, 2 any two values of variable- 

max maximum value of variable 


4 

cr 


x 


y 
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critical value of variable, value at which 
Mach number of 1 is reached at some point 
in flew 


partial 



partial 


derivative with respect to 

•N./ 


example : ^ ^ ~ 

derivative with resoect to 





THEORY 


The following is a condensation of the theory and mathematical 
development of the hydraulic analogy as given by Preiswerk (reference 4). 

Two assumptions are made in the mathematical development: 

(1) The flow is irrotat tonal. 

(2) The vertical accelerations at the free surface are negligible 
compared to the acceleration of gravity. The pressure in the fluid 

at any point therefore depends only on the height of the free surface 
above that point. 

The analogy between the flow of water with free surface and the 
flow of a compressible gas may be obtained by setting up the energy 
equations for each, From the energy equation for water, the 
velocity is 


and 

V 2 = 2g(d 0 - d) 

(la) 


^max = \/ 2gd 0 

(It) 

The corresponding 

equation for the velocity of a gas is 



V 2 = 2g(h 0 - h) = 2g Cp (T 0 - T) 

(2a) 


and 


^max ~ \Z 2 £ h o “ \/ 2 C c p- 0 


(2b) 
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Therefore, if the ratio v/v^^ for gas is equated to vA’max 
for water, 


d o " d _ lx o ~ h - __o 


- T 


ho 


or 


d _ T 

T’ ** t" 

a o O 


(3) 


The equation of continuity for water is 

u(ud) ->(vd) 

ir * = j 

The continuity equation for tvo -dimensional gas flov is 


S(up) d(vp) 


“ + TF 


= 0 


ox 


w 


(5) 


From equations (4) and (5) , a further condition for the analogy 
may “be derived 

d = P 
<k> p o 

Since for adiabatic isontropic flow in the gas 

1 

£_ - (X'7' 1 

Po '\ T o ! 


( 6 ) 


(7) 


and since, in the analogy, from equations (3) and (6) 


P 

p o 



T 


m 

-o 


m 

T 

x o 



then 
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and, therefore, the analogy requires that 7 
From the relation 

\7 


= 2 . 0 . 


p. = (p_Y 

p 0 ypoy 


it may he seen, since 7 = £.0, that 


P ./4f 

p 0 \ d °/ 


(8) 


The velocity potential for water is given hy the equation 


^oc 1 1 - ---■ J + i 


A 


\ fA 


77 


y 




- 2 <L~ = 0 

£&/ ' S& 


( 9 ) 


and the corresponding equation for a gas is 



Therefore, for identical expressions, 

o 

gd _ a7 
2gd c ‘ %h 0 


( 10 ) 


( 11 ) 


Freni equations ( 9 ) and (10), the velocity i/gd in the liquid 
flow is seen to correspond to the velocity of sound in gas flow. 

In a subsequent section entitled "Discussion," the value of ygd 
is shewn to be the velocity of propagation of surface waves, the 
wave lengths of which are large in comparison to the water depth. 

V 

The ratio of -- — in the liquid flow corresponds to the Mach 
\fed 

V . 

number — • m the gas flow. 

fct 

If the velocity of the liquid flow is less than y/g3 (M< l), 
the water is said to he_^' streaming." If the velocity of liquid 
flow is greater than ygd (M > l) , the water is said to be 
"shooting. " 
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In shooting vat or under certain conditions, the velocity ox" the 
flow may strongly decrease for short distances and the depth may- 
increase. An unsteady motion of this type is called a hydraulic 
jump. Hydraulic jumps of small intensity are propagated with 
the velocity \/gd. 

The analogy may he summarized as follows: 


Significant quantities 
and charact eristics 
of two-dimensional 
compressible gas flow, 
7 = 2 

Corresponding' values 
in analogous 
liquid flow 

T 

Temperature ratio, =j- 

1 o 


Water-depth ratio, •'x- 


Density ratio, fi-- 
80 


d 

Water-depth ratio, 

d o 


Pressure ratio, P— 
Po 

/? 

Square of water-depth ratio. 

( d '\ 2 

\ 4 o j 

Velocity of sound, a = 1 

Wave velocity , ^/gd 


w 

Mach number , - 
7 a 


Mach number, 

\ / cel 


Subsonic flow 


V t: 

Streaming water 


Supersonic flow 


Shooting water 


Shock wave 


Hydraulic jump 



The Mach numbers in the liquid flow may easily be computed. 
If equation (la) is substituted in the expression for stream 
Mach number 


’! S 


the expression becomes 


M 



1/2 


(12a) 


8 


MCA TN No. 1185 


The points for M = 1 are located vhere the depth is two-thirds of 
the total head as is easily verified from equation (12a) . The 
local Mach numbers are computed from 




U(do - d 1 ) 


1/2 



(12b) 


The pressure coefficient, a significant quantity in compressible 

/il'f 

gas flow, m ay "be computed for liquid flow. If ( — } is substituted 

Pq W 

for — in the expression for pressure coefficient 

P? 


the expression becomes 



P = 


V d s/ ‘ 1 


A oV 

W 


(13) 


For a compressible gas 


F = x 

c ntf 


a + 


X 

7 - 1 pY" 1 

M } 

2 s / 


- 1 


and with 7 = 2.0, as required by the hydraulic analogy, the 
compressibility factor is exactly 



When the analogy is applied to tho study of air flow, accurate 
quantitative results will not be obtained because, for strict 
agreement between water flow and gas flow, 7 must equal 2.0, 
whereas for air 7 is 1.4. The relation between pressure ratio end 
local Mach number for the two values of 7 is shown in figure 1. 

A given pressure ratio corresponds to a higher Mach number in air 
than in the fictitious gas with 7 = 2.0. Tile corresponding depth 
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ratios d/d Q shown in figure 1 arc equal to tlie square root of 

the pressure ratios for 7 = 2.0. With a given stream Mach number, 
the absolute value of the critical pressure coefficient (fig* 2) 
is lower in air than in the gas with 7 = 2.0. The differences 
are not very large, however, and the flow phenomena observed in 
the water flow should be qualitatively the same as those occurring 
in the two-dimensional compressible flow of air. 


APPARATUS 


The tests were conducted in the vertical, return-flow water 
channel which was designed and constructed in the Langley 8-foot 
high-speed tunnel. Figure 3 3hows the flow circuit and plan view 
of the free water surface and also the orifice locations of the 
water channel. A motor-driven propeller of 1-foot diameter forces 
the water into the expanding section and through a screen into a 
large quieting section. The fluid is accelerated in the convergent 
section and flows through the test section downstream to the propeller 
so that continuous flow is maintained. The test section could be 
replaced by the nozzles shown in figures 4(a) to 4(c). The shape 
of the channel test section is given for comparison in figure 4(d). 

The entire floor of the entrance cone and test section is horizontal 
except for the diffuser which has a 1° slope- for a static depth 
of 1.5 inches, the maximum volume flow is approximately 940 cubic 
inches per second and the power required is l/lO horsepower. 

Water depths at points in the test section were obtained with 
the survey equipment shovm in figures 5 and 6. A block on which 
was mounted a vertical shaft tipped with a fine wire probe (o. 030-inch 
diameter) could be moved both parallel and perpendicular to the 
channel axis. The probe was capab3.e of 2^ -inch vertical travel con- 
tinuously measurable by a standard micrometer to the nearest one-half 
thousandth of an inch. A small neon lamp was mounted in series with 
the insulated probe and the water to provide a positive sign of con- 
tact between the probe and water level. 

A similar micrometer unit, movable along a single fixed rail, 
was located upstream along the floor center line to measure total 
head. Water depths could also be measured by static orifices on 
the channel walls and along the floor. Each orifice was connected 
by a valve to a common sump and burette- By opening only one valve 
at a time, individual pressures could be obtained in the burette by 
means of a vernier height gage fitted with a hairline mirror sight. 

The locations of these orifices are shovm in figure 3* 
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Photographs end motion pictures of the free surface could he 
obtained by cameras mounted above the test section. A Strobolux 
and Strobotac were used as the source of illumination. A circular 
glass plate, fitted into the center of the floor of the test section, 
jjermitted illumination of the water surfaces from below. 


METHODS 


The two-dimensional, models were mounted at the- center of the 
glass plate, one end flush with the plate and the other end extending 
above the water surface. Orifices were drilled in the model a 
fixed distance above the glass plate and normal to the model surface. 

From the theory of the hydraulic analogy, it is readily seen 
that accurate measurements must be made of the depths in all 
measurable points in the channel. For each operating speed two 
readings were usually obtained for each data point, one hy the 
burette and the other by the survey carriage. The burette reading, 
by moans of the hairline 3 ight and vernier, could be taken as 
quickly or as often as necessary without affecting stream conditions; 
where as, because of the formation of capillary and standing waves, the 
surface proto reading was good only at tho instant the probe touched 
the water. This effect was not serious at local supersonic speeds as 
the standing waves were then unable to move upstream from the point 
of water contact. When the disturbance was removed by breaking 
contact with the water surface, original stream conditions were 
restored in about 10 seconds, at which time another reading could 
ho made. 


The micrometer unit movable along a single fixed rail measured 
the depth at the upstream positions and tho depths at the first 
pressure orifice on the channel center line. From continuity 
considerations, the total head 



where b is the breadth of the channel and subscript 1 refers to 
the station at the first orifice on the center line. The stream 
depth dg was measured sufficiently far ahead in the test section 
to be unaffected by the presence of the model. Th9 total head, 
however, could he measured directly by the upstream micrometer in 
the vicinity of the boundary wall at the upstream end (fig. 3) . The 
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upstream micrometer was calibrated against the tost -section survey 
probe before each run. /II depths were measured from the level 
floor of the channel. The net accuracy of fill readings was within 
TO. 002 inches. 

Orifice pressures were used to obtain the water depth on the 
surface of the models and at the test-section side walls. This 
typo of measurement is necessary inasmuch as the capillary rise 
of water on all bodies that break the water surface -invalidates the 
use of the survey probe. The probe is of great value. 1 however, in 
obtaining surveys at points whore the influence of capillarity is 
negligible, which in most cases is about l/4 inch from an exposed 
surface. Extreme care had to be taken in all pressure measurements 
to hold capillary rise constant in the burette, as very slight 
concentrations of d3rt on the glass weald change the meniscus shape. 
The burette was, therefore, cleaned at regular* intervals with 
chromic acid, and a 10 percent solution of aerosol was placed cn 
the meniscus to reduce surface tension and capillarity. This 
procedure enabled consistent readings within the desired accuracy . 

The shadowgraph still photographs wero obtained by a r j- by 7 -inch 
view camera. The test setup and schematic diagram ace shown in 
figures 6 and 7(a) . In the initial trials it was found that the 
most convenient light source, with existing collimator lenses, was 
a E br obolux -Strobotac unit) ani b,> adjusting the flash frequency, 
various light intensities could be obtained. The refraction 
patterns formed by a vortex, hydraulic jump, and capillary wave 
are shown in figures 7(b) to 7(d), and from these patterns or 
combinations, a qualitative analysis of shadowgraphs may bo made. 

In order to reduce excessive pattern ref ract i cn, the ground-glass 
screen was placed as close to the water surface as possible. Motion 
pictures of wake frequencies and alternate trailing vortices were 
obtained by the same methods with no great difficulty . 


RESULTS 


In order to ascertain the nature and the applicability to air 
flow of data obtained in the water channel, tests wore mad© with flow- 
through nozzles and about circular cylinders. These tests wore- run 
at various depths and Mach numbers end with models of various sizes 
in the channel. The results are first presented without discussion. 

Plots of the water depths and water-depth ratios along the center 
line of convergent -divergent nozzles as measured by the test -section 
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mi creme ter unit aro shown in figures 8 to 10. The shapes of these 
nozzles are shown in figures 4(a) ^ 4(h), and 4(c), respectively. 

The curves labeled "Reynolds flow'' are computed froEi one -dimensional 
theory) that is, the velocities at all points in a cross section 
are assumed the sane in magnitude and parallel to the center line 
of the channel. The depth at any point is computed from the 
equation of continuity 


Y x Area = Constant 


(15) 


where 


Area = Breadth of the channel X Water depth 

The water-depth ratios are seen to increase with the increase in 
depth in figure 8. The nozzle in figure 10 was modified to have 
a very gradual approach section and the v water-depth ratios for 
this nozzle ere the same at any depth. 

The variation with total head of the water depths and water-depth 
ratj.es about a circular cylinder -as measured by means of the 
pressure orifices and burettes is shown in figure 11- In this 
case ths depth ratios decrease with increase in water depth. A 
comparison between micrometer and burette measurements is shown 
for. a C-inch nozzle in figure 12, where it is seen that the 
micrometer readings are greater than the burette readings and the 
difference increases with depth. 

Figures 13 and 1’:- show the variation with water depth of the 

shadowgraphs of the flow about a 2~-inch diameter circular cylinder 

at Mach numbers of 0.40 and 0 - 60, respectively. The corresponding 
variation of pressure at d = 90 0 is also shown. Additional 
shadowgraphs of the flow about circular cylinders of various sizes 
and with various water depths are shown in figures 19 to 19- The 

pressure distributions about the 2-f-incli circular cylinder in the 

o 

water charnel with a static depth of ^ inch are compared for various 

Mach numbers in figure 20 with the pressure distributions about 

a y-inch circular cylinder in air. The data for the circular cylinder 

in air were available from tests made in the Langley rectangular 
high-speed tunnel. 

With various sizes of circular cylinders in the test section the 
variation of the local Mach numbers along the channel wall for different 
stream Mach numbers is shown in figures 21 to 2b) and the corresponding 
stream depths d E and total heads d Q are given in figure 25 . 
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DISCUSSION 


Various factors influence the results obtained in the water 
channel and the effect of these factors must he ascertained 
before the channel can be, used for the investigation of the 
analogous gas flows. 

From an examination of figures 15 to 19 it is seen that 
concentric waves appear in front of the model. The fact that 
these waves do not appear in the schlieren photographs (fig. 26) 
of air flow taken at the Langley rectangular 1 high-speed tunnel suggest?? 
that these waves are not part of the analogy. 

A disturbance at some point in a liquid generally will give 
rise to two types of waves: short surface-tension waves called 

capillary waves and considerably longer gravity waves (reference 
5 and pp. 353-402 of reference 6) . The expression for the 
velocity of propagation of water- waves is 



(16) 


Figure £7 is a plot of IJ against X for water depths of 0.5 inch 
and 2.0 inches. This figure shows that waves cannot exist with a 
velocity of propagation less than approximately 0.75 foot per second. 
Disturbances of wave length less than that corresponding eo the 
minimum wave velocity ere termed capillary waves since they depend 
primarily on the surface tension of the fluid. Propagation velocities 
greater than this minimum velocity correspond to shorter capillary 
waves and longer gravity waves. 


If a = 0 and d « X, equation (16) becomes 

U = y'gd (17) 

and this velocity - the velocity of propagation of long gravity waves 
is the basic surface wave velocity. The other waves defined by 
equation (16) are purely water surface waves and are not considered 
in the hydraulic analogy. Figure £7 shows that changes in the 
water depth have little effect on the velocity of propagation of 
waves for wave lengths less than 0.1 foot. 

The wave lengths of some of the standing waves appearing in 
front of the models in figures 15 to 19 were measured from the 
photographs and were plotted in the inset of figure 27 against 
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the local velocity (equal, for standing waves, to the velocity of 
propagation) at the point where the wave occurred. This plot 
shoved that these waves were capillary waves. The check marks 
on the inset show the number of readings that fall on the 
plotted points. These capillary waves, although they have no 
part m the hydraulic analogy, have two adverse effects on the 
application of the analogy; namely, complication of the flow 
photographs and decrease in the accuracy of the depth measurements 
taken with the test -section probe. 

The results obtained in the water channel depend to a large 
extent on the depth at which the tests ere run. The effect of 
depth on the formation of standing capillary waves is shown in 
figure 28. If high stream Mach numbers without capillary waves 
in the stream are to be obtained, a shallow depth is prescribed.. 

The effect of depth on the flow patterns and pressure 
coefficients about circular cylinders is very pronounced (figs. 13 
end lh) . The pressure coefficient rises to a maximum value 
ac dgt 2 0.9 inch. The decrease in pressure coefficient 
at lower depths is believed to be due to the affect of bottom 
boundary layer; the decrease at higher depths is duo to the effect 
of the vertical accelerations. Comparison of figures 13 ana lb 
with figure 26 indicates that the water flow most nearly corresponds 
to air flow at depths that arc between 0.79 inch and 1.0 inch. 

The water is accelerated from the total head where the fluid 
velocity is zero to a value of zero acceleration and maximum 
velocity in the test section- The motion of the water through 
this cycle is actuated by the force of gravity which is thus the 
important factor in the hydraulic analogy. If an obstruction is 
placed in the path cf a fluid at constant velocity, local accelera- 
tions of flow must take place in the field about the obstacle. If 
these accelerations arc large, the vertical components, which have 
no part in the analogy, are no longer negligible in comparison with 
the acceleration of gravity; and a distortion cf the flow results. 

This distortion in such that in a region of deceleration such as 
the region in which the flow is approaching: stagnation, the depth 
is less than that required in the analogy; and the depth Indicated 
by the pressure in the static orifices is greater than the true 
depth. In a region of acceleration the opposite effects occur. 
Differentiation of equation (la) shows that the acceleration 
(in steady flow) is proportional to the gradient of the depth. 

The slope of the free surface is therefore a moaouro of the 
accelerations . 

If, for two cases of the same stream Mach number the water 
depth is increased, the slope of the free surface is likewise greater as 
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shown in the following table in which for chosen values of M end 
d Q the value of d s is computed from equation (12a) : 


M 

d o 

d s 

d c/ d o 

1 

0 

V(in./sec ) 

0.926 

.926 

1.000 

1.100 

0.700 

.770 

0 0 
00 
ir— 

c 

0.300 

.33c 

10.78 

11.29 


The term (d c - d s ) is an indication of the slope of the free 
surface; the higher value corresponds to a steeper slope and thus 
to greater vertical accelerations- An illustration of the fore- 
going statements is given by the results of the nozzle tests in 
which; as the water depths decrease, the depth ratios approach 
those calculated by Reynolds flow (figs- 8 and 9); all tests were 
made at choking Mach number. Binnic and Hooker observed the same 
effect (reference 3) • If the slope of the free surface is small 
and if the nozzle is so shaped that the acceleration in the entrance 
can be kept sufficiently small, no variation of water-depth ratios 
occurs with changing depth (fig. 10) • The same effect is shown in 
figures 11 and 12 as measured by static orifices find burettes. 

The effect of vortical accelerations on depth ratios determined 
by means of static orifices is opposite to that obtained by the 
probe. It may be assumed therefore that this adverse effect would 
bo negligible when the two methods of reading pressures coincide. The 
readings taken at the minimum section of the 2-inch nozzle (fig. 12) 
when the depth was reduced to 1.3 inches coincided. If the free 
surface slope had been more abrupt (smaller model or model having 
a sharp pressure rise) , a still shallower depth would have been 
required to attain this condition; therefore large models and small 
depths should be employed. The minimum depth is limited, however, 
by the effects of tho boundary layer on the floor of the channel. 

The development of tho boundary layer produces a velocity gradient 
along the channel similar to the gradient in a wind tunnel. This 
effect can be compensated in the same way by diverging the walls or' 
in the case of the water channel by sloping the floor downward in the 
direction of flow. An effect for which compensation does not appear 
feasible is the thickening and thinning of the boundary layer in 
regions of decelerating and accelerating flow about a model, behavior 
which causes a distortion of the velocity distribution such as would 
correspond qualitatively to a thinning of the model. These boundary 
layer effects are particularly serious in the water channel because 
of the appreciable boundary layer displacement thickness (found to 
be of the order of 0.1 in.) in comparison with tho total depth of 
the fluid. The boundary layer effects could obviously be minimized 
by increasing the water depths, both because of tho increased ratio 
of total depth to boundary layer displacement thickness and 
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"because of the increased Reynolds numbers corresponding to the 
highor velocities with given Mach lumbers . Thi3 requirement is, 
unfortunately, inconsistent with the minimizing of the effects 
of the capillary waves end of the vertical accelerations. 

The analogy is further limited by the fact that it applies for 
a hypothetical gas having a value of 7 = 2.0; the effect of this 
limitation has not been completely determined. The effects shown 
in figures 1 and 2 are not very large, however, and can be taken 
into account in the interpretation of" data taken in the water 
channel. Other investigations suggest that the influence of the 
value of 7 in subsonic compressible flow is not great . Kaplan 
in reference 7 thus found that to the third approximation the 
effect of the value of 7 in the c ompr e s s i bl 0 flow up to the 
critical speed was nevlisible and von ICarman, reference 8, obtained 
a widely used expression for change of pressure coefficient with 
Mach number by using the assumption 7 = -1. For supersonic flow, 
the characteristics curves are considerably influenced by the value 
of 7 . The ratio of maximum velocity to the velocity of sound is 
greater with 7 = l.k than with 7 = 2.0, and for a given change in 
flow angle, the velocity change, as determined by the characteristics 
method, is also greater. 

The Reynolds numbers in the water channel are calculated for 
complete submersion and are very low; the maximum value of E for 

the 2^-inch-diameter cylinder with a static water depth of 0.75 inch 

is 15,000. The comparison of the water-channel and wind-tunnel data 
(figs. 15 to 19 and 26 ) showed that, although the cylinder tests 
were made at considerably higher Reynolds numbers in air than in 
water, the same type of flow occurred in both cases; the main 
consideration is agreement of the flow patterns rather than 
identical Reynolds numbers. 

In a wind tunnel, a. choking or maximum stream Mach number occurs 
when local Mach numbers of 1.00 extend across the section between 
the model and tunnel walls. A similar effect was noted in the 
water channel. Figure 25 shows that the stream depth decreases and 
the total head increases with Mach number up to a maximum value 
that depends on the size of the model tested. At this choking 
condition, both the stream depth and the total head Increase with 
an additional power input. Figure 29 shows the maximum Mach number 
for various ratios of cylinder diameter to channel width. Variation 
of the data from tho theoretical curves might be expected inasmuch 
as the theoretical values wore computed by one -dimensional theory 
(Reynolds flow); whereas, tho actual flow is two-dimensional. 

Maximum Mach numbers that occur in wind tunnels, however, agree much 
more closely with the one -dimensioned theory than do those shown 
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in figure 29 (see reference 9) , and it therefore scens likely" that 
the divergence is due to other factors, such as "boundary -layer effects 
and vertical accelerations, rather than to the cvo-dimensional nature 
of the flev. 


The disturbance along; the channel wall with various size 
cylinders in the test section is shown in figures 21 to 2k. The 
wall disturbance is an indication of the severity at high Mach 
numbers of channel- wall or wind-tunnel -vail interference. As the 
stream Mach number approaches 1.00, the wall disturbance increases 
until at the choking condition the interference becomes large even 
for the smallest cylinder. For the 1— -inch- diameter cylinder 


(fig. 22), the point of sonic velocity at the wall is located 
approximately 3 diameters behind the cylinder . This result is in 
agreement with the flow photographs (figs. 15 and IT), which show 
a gull -wing-shaped shock starting from the region behind the 
cylinder and extending downstream toward the walls. Figure 30 
shows the lines of constant Mach number in the flow field about 
a 6 -inch-diameter cylinder at choking Mach number and shows clearly 
the supersonic flow occurring behind the cylinder. 


A comparison of the shadowgraphs (figs. 15 to 19) with schlieren 
photographs of flow about circular cylinders in air (fig. 26) shows 
that the two flows are very similar. The typo of flow is 
the same in the' two cases and is the type that is characteristic 
of Reynolds numbers somewhat below the critical value. At Reynolds 
numbers between 50 and 350,000, an unstable condition is set up 
behind a circular cylinder in which vox'tices are shed alternately 
from each side. Laminar separation occurs at a point approximately 
80° from the forward stagnation point and a vortex sheet is formed 
which extends downstream and finally rolls up into a large vortex. 

The voruices so shed arrange themselves into a Xarman street. (See 
reference 6, pp. 21T-218.) Figures 16(a), 13(a), 19(b), and 26(a) 
show this type of flow, which is the same in the water channel as 
in the air flow. Because of the occurrence of separation, the actual 
pressure distribution is different from that calculated by potential - 
flow tlieory . The actual pressures over the forward part of the 
cylinder are higher and those over the rear are lower than those 
calculated by theory. The negative pressure peak is greatly reduced. 
Figure 20 shows the pressure distributions about circular cylinders 
in air and in the water channel. The quantitative results obtained 
in the wate v> channel are very close to those obtained in air 
(figs. 20(a) to 20(c)). This close agreement is believed, however, to 
be largoly fortuitous, resulting from an interaction between the bottom 
boundary layer effects and the rather large channel -wall interference 
in the water channel. At the choking Mach number M = 0-T1 in "the 
water channel (see fig. 20 (d)) the - interference effects exceed the 
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relieving effects of the thinning of the “boundary layer and the 
negative pressure coefficients near 8 *= 90° exceed those obtained 
in the wind tunnel* 


Because of the unsteady nature of the flow , instantaneous 
velocities greater than those indicated "by the pressure distributions 
occur near the surface of a cylinder. The critical speed may 
therefore be expected to be lower than indicated by the pressure 
distributions, an effect that has been observed in air * Figure 2 
shews that with a given peak negative pressure coefficient, the 
critical Mach number should be somewhat lower in the water channel 
(7 = 2.o) than in air (7 = 1*4)* 

Another effect of the unsteady flow is the asymmetrical 
development of shock* When the flow closes in on one side of the 
cylinder, the induced velocities increase and may exceed a Mach 
number of 1*00* A shock wave may thus form on one side while none 
exists on the other, or a more intense shod: may exist on one side 
than on the other. Such asymmetrical shock patterns are shewn in 
figures 16(b) to 16(d), 17(b), 18(b) to 18(d), and 26(b) to 26(d)* 
Exactly the same phenomenon occurs in the water channel as in the 
air flow* In seme cases, the shedding of these waves alternately 
from each side of a cylinder has been observed in the water channel 
when the stream Mach number was only a little above the critical value* 
Such behavior is to be expected if on one side of the cylinder the 
velocities first exceed and then fall belcw the speed of sound as 
the flow closes in and then breaks away from the surface. With 
increase in Mach number, a strong disturbance originates at the 
edge of the wake approximately 1 diameter behind the cylinder and 
extends into the field of flow (figs* 15(d), lo(d), and 26(d)). This 
disturbance oscillates with the wake, which is still unstable, and 
in its incipient stages alternates from one side of the model to the 
other. At still higher speed, the flow closes in behind the cylinder 
so that the cylinder has the appearance of a streamline body with a 
strong gull -wing- shaped disturbance at its trailing edge (figs* 15(e), 
15(f), 16(f), 17(f), 18(f), 19(f), 26(e), and 26(f)). Thesd 
features appear to be essentially the same in the water channel as 
in the air flow* 


EQUIPMENT AND FUTURE DEVELOPMENT 


Measurement of depth at the surface of models and at fluid 
boundaries are best made with pressure orifices mid burettes 
because of the capillary rise, whereas the micrometer method is 
suited for measurement of field depths; both methods agree very- 
well provided the vertical accelerations are small* The burette 
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method is best for obtaining average pressures in -unstable flow. A 
method of photography has been developed that is simple and completely 
satisfactory. The shadowgraphs of the water flow are strikingly 
similar to the schlieren photographs of air flow. 

The value of the water channel lies mainly in the low cost 
and convenience of operation, field surveys are simply made, and 
various features of the flow such as turbulence, vortices, separation, 
fmd shock formation are easily observed and photographed. Streamlines 
about a model are easily obtained by inserting streamers of dye in 
the water ahead of the model. 


A larger channel especially designed to minimize boundary layer 
effects and to secure uniform flow in the test section is desirable. 
The channel should be large enough to permit tests at Reynolds 
numbers above the critical value. A larger channel would also 
be advantageous in reducing the adverse effects of the boundary 
layer or vertical accelerations, or both. 

Additional investigations, both theoretical and experimental, 
are needed in order to determine the corrections necessary to 
convert quantities obtained from the water flow to the values 
characteristic of the air flow about corresponding bodies. 


CONCLUSIONS 


An experimental apparatus end technique have boon developed 
for the investigation of the analogy betwoon water flow with a 
free surface and two-dimensional compressible gas flow (hydraulic 
analogy) ; a preliminary investigation has been made and the results 
of an application of the analogy have been presented. The .following 
conclusions are indicated from this work: 


1. The hydraulic analogy provides a very inexpensive and 
convenient means of investigating high-speed two-dimensional air 
flow. Tho flow may be observed and photographed, and surface and 
field measurements may be easily obtained. Reasonably satisfactory 
agreement was found between the water flow and air flow about 
corresponding bodies, although considerable work in both theory 
and experiment is needed in order to convert with quantitative 
accuracy from the •'.•miter flow to the flow in air. 
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2. With a larger channel, difficulties due to vertical 
accelerations and subcritical Reynolds numbers might he overcome. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , August 21, 1946 
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Figure 3. ~ Schematic v/err of rvafer channel rv/th seef/ons. 
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Fig. 5 



Figure 5.- Test section and survey equipment seen from downstream end of channel. 


j 



NACA TN No. 1185 


Fig. 6 



Figure Approach and test section showing equipment for photography. 
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Fig. 16 
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